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FOREWORD 


This technical report concludes Contract No. NAS 9-14349, "Zero- 
G Flight Test of a Gauging System, " and covers the experimental flight 
testing of a nucleonic quantity gauging system conducted on-board a KC- 
135 Zero-G Aircraft, The report was prepared by Frank E, Bupp of 
TRW Defense and Space Systems Group, One Space Park, Redondo Beach, 
California. The author is grateful to numerous personnel for assistance 
in completion of this program, particularly to Mr. M. McFarlin of TRW, 
who was responsible for the development of the flight system software 
and data reduction and analysis, and Mr. R. Cardon of TRW for electron- 
ic support. The author also expresses his grateful appreciation to Mr. 

J. Alexander of the NASA/ JSC and to Mr. D. Griggs, Chief, Zero-G 
Project Office, and the members of the flight crew for their contributions 
and dedicated participation in the conduct of the flight tests. 

This report contains no classified information extracted from other 
* classified reports; this report was submitted in January 1976. 
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SECTION I 


INTRODUCTION 


Since November of 1970, a nucleonic gauging system for gauging 
propellant in tankage under zero-gravity conditions has been under devel- 
opment at TRW Systems. This gauging system is based upon the principle 
of "X-raying" the interior of the tank wherein the absorption of the diag- 
nostic radiation by any intervening propellant is utilized to determine the 
propellant mass. The technological feasibility of this diagnostic technique 
was first demonstrated successfully in June of 1970 in the form of a bread- 
board model installed on battleship tankage under one-g conditions. In ihe 
interest of minimizing cost, this early feasibility model was interfaced 
with a large, general purpose computer for requisite data processing. The 
results obtained with this early breadboard model were, to say the least, 
promising, particularly in the high levels of accuracy that were achieved. 
Consequently, a complete flight-type system including an integral micro- 
computer was next designed, fabricated, and subsequently tested in 
November of 1973. This fli fe ht-type system, described in detail in Ref- 
erence 1, was highly successful and achieved total system accuracies of 
0. 3 5 percent and better; of fundamental importance to this high level of 
accuracy, however, was the system software (i. e. , requisite to the mass 
computation). This report, therefore, documents the nucleonic gauging 
system software for the computation of propellant mass; this software 
a Zero-G Gauging System, 11 and subsequently modified to meet the require- 
ments of zero-g flight testing under Contract No. NAS 9 - 14349 . 


1-1 



REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


SECTION II 

GENERAL DESCRIPTION OF THE 
NUCLEONIC GAUGE CONCEPT 


In its most general form, a nucleonic gauge is comprised of an array 
of gamma ray emitting radioisotope sources positioned on one side of a 
(propellant) tank, and an array of radiation detectors positioned on the 
other side of the tank opposite the sources; the concept is illustrated in 
Figure 2-1. Although all the elements and subsystems of the gauge are 
external to the tank, gamma rays that are emitted toward the detectors 
will penetrate the tank walls and be partially absorbed by the intervening 
fluid (fuel) within the tank. The amount of radiation which emerges on 
the detector side of the tank is then proportional to the mass of the inter- 
vening propellant. 


In practice, there are two basically different methods of mechaniz- 
ing a nucleonic gauge: one mechanization is based on the natural — i. e. , 
exponential — gamma ray absorption characteristics of the propellant 
while the other is based on obtaining a linear gamma ray absorption rela- 
tionship, The latter technique is accomplished by selectively employing 
a window discriminator in the detection circuitry to accept only those 
photons which fall within a narrow energy range, thus making the gamma 
ray absorption a linear function of the absorber thickness. Although the 
linear system is the easier of the two systems to construct (at least from 
a hardware point of view), it is severely limited to tank dimensions the 
order of 8 feet or less; this limitation precluded the linear system from 
further consideration for the envisioned application and hence it was not 

subjected to rigorous study during 
this program. This report dis- 
cusses only those systems based 
on exponential absorption. 

For a propellant gauge 
employing the exponential principle 
of gamma ray absorption to operate 
properly under low Bond number 
regimes and/or zero-g conditions, 
it must first be insensitive to the 
fuel orientation within the tank. 

This is accomplished by judiciously 
distributing the detectors and 
radioisotope sources opposite each 
other and assuring that the detec- 
tors and sources are configured 
into a suitable array which covers 
a substantial portion of the entire 
projected area of the tank. The 
source- detector array is absolutely 
necessary unless the fuel is con- 
strained by some means (such as 



Configuration 


2-1 




‘ * . 



I 


I 


I 

! 


centrifugal force) into a known geometric distribution. For example, it 
would not be possible to use a single radioisotope source in the (central) 
interior of the tank and an array of radiation detectors on the exterior. 

If such a system were used to gauge propellant that may assume a random 
distribution (as in zero-g), a large error in propellant quantity measure- 
ment would result as is clearly envisioned in one case where the propellant 
is distributed along the walls of the tank, whereas in another case it is 
concentrated around the source in the .enter of the tank. The mean pro- 
pellant thickness is the same in both cases. S’ nee gamma ray absorption 
is a function of the thickness of the absorber, there will be equal amounts 
of gamma radiation absorbed within the propellant in both cases, and 
equal amounts of radiation will reach the detectors despite the fact that 
the first case has perhaps nearly an order of magnitude more propellant 
than the latter case. Needless to say, an attempt to use such a system 
for gauging the propellant quantity will produce a very large error. 

To overcome this problem of random fuel orientation, the sources 
and the detectors are configured into arrays and distributed on a (nearly) 
planar surface, thereby creating a uniform gamma ray field within the 
boundary of the tank. In this gamma ray field, an element of propellant 
may be moved anywhere within the boundary of the tank and, regardless 
of its position, will absorb the same amount of gamma radiation. 

Due to limitations on the number (and distribution) of sources and 
detectors that can practicably be employed, nucleonic gauging invariably 
depends upon discrete sampling of the tank interior; a typical configuration 
is illustrated in Figure 2-1. 


As shown, the gauge consists of individual gamma ray radioisotope 
sources, each of which is external to the tank and collimated to produce 
radiation in a narrow cone. On the opposite side and also external to the 
tank, a single detector is positioned in the center of each of the cones of 
radiation. Each source is sufficiently collimated so that it essentially 
illuminates only its opposing detector. The combination of each source 
and its opposing detector comprise a "source-detector pair", thereby 
providing a discrete "sample column" through the interior of the tank. 

By configuring a number of such source- detector pairs into an appropriate 
"sampling array", the entire contents of the tank can be sampled. 

In general the requirements for a radioisotope source for propellant 
measurement are high energy penetration capability to permit source - 
detector installation external to the tank, a comparatively long half-life 
to minimize calibration, and a source strength consistent with the accu- 
racy and resolution requirements without compromise to personal safety. 


2. 1 PRINCIPLE OF OPERATION 


.th 


In the nucleonic gauging system, the mass sampled by the i— sam- 
ple column due to the exponential attenuation of primary radiation from a 
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radioisotope source when an absorber and/or fuel intervenes the gamraa^ 
ray beam will be proportional to the number of photons received at the i — 
detector, , according to: 




(2- j. i 


where 


N 


o 


number of counts received in the absence of an intervening 
mas 8 


-p- = mass absorption coefficient 

P = density of the intervening mass 

Z = thickness of the intervening mass 

In the event of a tank containing both fuel (liquid) and gas (ullage), 
the number of counts received by the i& detector is simply 



( 2 - 2 ) 


where the subscripts L and G refer to the liquid and gas, respectively, 
within the ith sample column. 

The mass, M., in the ith column is given by 

M i - ,W i )A i (VV G Z G f ) 12 
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where is the area associated with the i& sample column and W^ is the 
weighting factor associated with the itk source- detector pair. In the 
typical situation where the sample columns are parallel 


A. * 

1 


and W. * 1 for all i 
1 


(2-5) 


where Ax is the total projected cross-sectional area of the tank and k is 
the number of sampling columns. 

Thus, the total mass sampled by the k sample columns is 


M 



M. = 




( 2 - 6 ) 


or: 



(2-7) 
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Figure 2 -2. Two-Dimensional Area- 
Weighting for Comput- 
ing Propellant Quantity 


This technique for com- 
puting propellant quantity is 
graphically illustrated (in 
two dimensions) in Figure 2-2. 

It is now noted that the 
number of sample columns are 
finite and, hence, tl e "cal- 
culated" mass will not (nec- 
essarily) be precisely the 
' actual" mass. Also, the 
calculated mass will depend 
upon the configuration (spatial 
distribution) of the sample 
columns and the tank size ami 
geometry. What is required, 
then, is a conversion — or, 
more correctly — a calibration 
curve from which the actual 
mass can be determined from 
the calculated mass, i. e. , 


TUTR' ‘f'TJOIBILnY OF THE 

■V < \ V RACE IS POOR 


2-4 





In practice, Equation (2-8) is solved with a computer simulation 
wherein the propellant can assume random orientations within a propel- 
lant tank under zero-g conditions. (See Reference 1. ) However, Equa- 
tion (2-7) is very easily mechanized with a small capacity digital computer 
and, a priori, so too are the functions represented by Equation (2-8) for 
one-g conditions. 



SECTION r. 


GAUGING SOl’TWARE DESCRIPTION 


A major consideration for any flight article is the minimization of 
both size and weight. Therefore, considerations for a system controller 
computer were in the field of mini- or micro- computers. A Control Data 
Corporation, model 469 mini- computer was chosen because of its high 
reliability in applications wi.ere small size, light weight, and low power 
requirements are a must. (A complete description of the computer can 
be found in Section 6. ) 

The software used for the gauging system is referred to as the 469 
assembly language. This language is similar in application to any assembly 
language but unique to the 469 computer. A ccmplete explanation of the 
language is presented in References 2, 3 and 4. 

Figure 3-1 presents a flow diagram of the overall program structure. 
The program can easily be separated into three operating modes, each 
independent of the other except for data linkage. These three modes are: 
Program Initialization, Mass Calculation, and Data Verification. 

3. 1 PROGRAM INITIALIZATION 

The -'omputer used for this system is a CDC 469 mini- computer 
with 20, COOg words of memory. Of this memory, the first 2000g words 
are defined by hardware as the DRO (destructive readout) region commonly 
referred to as the scratch pad. Of this 2000g words of scratch, the first 
4008 words are referred to as page zero, and the first 20g locations of 
page zero are the processor registers. Refer to Figure 3-- for a graphic- 
description of the memory allocations. 

The DRO region has three characteristics that play an important 
part in the development of the program. First, this region may be 
destroyed randomly at power-up and power-down or with the depressing 
of the master clear button. This makes it impractical for program stor- 
age. Since hardware defines a starting location of 400g, the first function 
of the software is to transfer control to the NDRO region (20008), which 
is preserved during these critical times. Secondly, the DRO region can 
be written into during program execution. (The NDRO region cannot be 
written into without first depressing the NDRO button on the programmer's 
console. This prevents inadvertently writing over the original program 
during execution. ) Therefore, the DRO region is used for temporary 
data storage. Thirdly, page zero has the unique ability to be addressed 
directly for read and write instructions; other pages must be referred to 
by indexed or indirect addressing. This ability to read or write with direct 
addressing saves both execution time and program storage. Page zero 
is best used for storage of any constants that are used frequently. 
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Figure 3-1, Program Startup 


The first function of the software is to transfer control to 20008 
(the start of the NDRO region); the second function is to zero the scratch 
pad from one location 20g to 17778* All the frequently used constants 
can now be transferred to page zero. 

As mentioned previously, hardware defines a starting address of 
4008- At power-up, the P register (program counter) will therefore 
display 4008. Initially, control is transferred to START (20008) by an 
unconditional jump (UJ), which is stored at core locations 4008 and 401 a. 
Since this instruction is in the DRO region, it may be inadvertently 
destroyed; it is rewritten to lower core from the write protected region 
at each program startup. 
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Certain options have been incorporated into the software to operate 
various peripherals. If the data are to be retained for further analysis 
via either mag-tape or teletype, the corresponding flag must be set accord- 
ing to Table 3-1. The system is constructed to run completely independent 
of any peripherals, or with the mag tape and/or teletype connected. One 
other option is the use of the Signal Conditioning Electronics for data input 
or to use constant data. Using constant data is useful for program check- 
out, though once the program is functioning properly, the inputs from the 
SCEs should always be used. 


Table 3-1. System Flags 



Core 



Variable 

Location 

Value 

Description 

GOMT 

2023 g 

0* 

Perform mag -tape write 


1 

Skip mag - tape write 

GOTEL 

2025 g 

o* 

Perform teletype write 


1 

Skip teletype write 

KFLAG 

2027g 

0* 

Use SCE input 


1 

Use constant inputs 


* Nominal values 


If indeed, the mag-tape is to be used, one other function is to be 
performed by the initialization portion of the software. The mag -tape is 
positioned at the beginning of tape and a dummy file is generated. This 
dummy file is used by the search and read routines for the mag tape. 

Control is now transferred to the teletype, and the flag "Ready" is 
displayed. At this time any descriptive information may be entered onto 
the mag tape. Usually at this time, the date, fuel loading, and any partic- 
ulars of the run are entered. Input is terminated with the $ key, and con- 
trol is again transferred to the console. The program is now in a pause, 
waiting for a level three interrupt issued by hardware to acknowledge the 
beginning of the mass calculation sequence of the code. 

3.2 MASS CALCULATIONS 

This portion of the program will continue uninterrupted until such 
time as an input from the console is received. The total propellant mass 
is basically generated by a polynominal curvefit. The input from the 
SCE's records the actual number of counts received by the detectors for 
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any 1- second interval. These count data can then be directly converted 
to a ray travel distance through the intervening mass by the following 
equation: 



where: 

= ray travel distance through the intervening mass (cm) 
fi = propellant absorptivity (cm~*) 

N. = counts received via channel i 

i 

= counts received for an empty tank via channel i. 
i 

The Z’s can then be used to calculate the propellant mass within the 
tank by means of a curve-fitting subroutine stored in core in the 469. This 
curve -fitting subroutine is referred to as the calibration curve(s) and is 
previously determined off-line in a large scale scientific computer. Es- 
sentially, a Monte Carlo simulation of the propellant behavior in zero-g 
conditions for the given tank/ gauge get metry is employed to rtlate the 
(average) Z’s to the actual mass within the tank. This Z-mass relation 
ship is (invariably) in the form of a polynomial, and it is the degree of 
this polynomial and its coefficients which are stored in the 469 core. 

Note that this polynomial (the calibration curve) is determined uniquely 
for each particular tank/propellant/ gauge -configuration. 

For each source/detector pair there is both a Z and a mass calcu- 
lation made based on its respective calibration curve. After the masses 
have been computed, a failure analysis must be performed to check their 
validity. The four masses (there are four source/detector pairs -channels) 
are compared against one another, and the differences are recorded. If 
one mass consistently has a large difference when compared to the others, 
that channel is considered inoperable, and the respective mass is weight- 
ed out of the total mass calculation. The system can operate with the loss 
of one channel, but the loss of two or more channels wall result in the 
shutdown of the entire system. 

In the event of system shutdown the display is blanked, but the sys- 
tem will continue to operate. Inputs are still accepted and masses gener- 
ated, and the failure analysis still performed. If for any iteration the 
masses have stabilized sufficiently to pass the failure analysis, the display 
will monitor this new mass, and the system will begin to operate again. 

The mass calculations will continue indefinitely until instruction from the 
console to halt the mass calculation sequence is received. 
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3.3 DATA VERIFICATION 


After each iteration of the mass calculations, the console is checked 
for an input. At the sign of the first non- zero input the system will dis- 
continue its mass calculations, write an end-of-file mark onto the mag- 
tape, and pause for further instructions from the console. At this point 
the operator has two choices: first, he may push the STEP, MASTER 
CLEAR, and RUN buttons, in that sequence, to start the program from 
the very beginning, or he may verify that data have been entered and writ- 
ten to the r: ag-tape. By changing the input from the data register of the 
console, the program goes into the data verification portion of the soft- 
ware. By data verification, the system verifies that data have been taken, 
not that the data are correct. To do this, the mag-tape is rewound one 
file, and the data acquired during the last phase are displayed on the tele- 
type. Six columns of information are displayed, all in octal. The first 
four are the count data accepted from the four source -detector pairs. 

This information will enable the operator to see if the four channels are 
operating correctly. The tifth column will be the total mass of the system 
generated by the zero-g gauging system. The sixth column will be the to- 
tal mass generated by 4 'he Rho-V gauging system (not implemented at this 
time). Cnee again these quantities are displayed in octal and must be con- 
verted to decimal for comparison with the mass displayed by the system 
hardware. 

After completion of the data dump to the teletype the system will 
restart itself by displaying "ready" on the teletype. (Note: the system 
should be allowed to display all the data to insure correct system operation. 
Early termination will result in the loss of data. ) At this point, the 
program has been re- initialized, and the operator can enter new run infor- 
mation and start the mass calculations again. 

3. 4 CHANNEL INPUT DATA FORMAT 

The data read from the accumulators must be made available to the 
469 computer in a very specific manner. The total counts received by 
each individual accumulator may range from 100 counts (a full tank) to 
a maximum of 2 million counts (an empty tank). Since a normal, single- 
precision data word represents 15 bits of information, a maximum data 
value of 2*6 _i (65535) implies double precision inputs are necessary. 
Therefore, two words of information, and thus two input channels, must 
be used for each accumulator. 

Double precision is correct only in the sense that two data words are 
used for input. In reality, a modified word organization is utilized. The 
entire second word is not needed, and status information is inputted as an 
additional 4-bit piece of information as nhown in Figure 3-3. Notice that 
the "sign" bit in the second word is used as a data- bit along with the next 
adjacent five bits. This enables a maximum reading of 2 21 -l (2, 097, 151) 
counts to be processed. The next four bits represent a status flag, indicat- 
ing the operational status of the four external accumulator devices. 

Modified program execution must be made as a function of the failure of 
any or all of the accumulators. The last six bits are set to zero and serve 
no particular function at this time. 
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DATA INFORMATION 


T 


SIGN 

BIT 

+ 


WORD 1 

(EVEN INPUT CHANNELS) 


CHANNEL 


2 5 ? 4 2 3 2 2 2* 2° • 4 3 2 I • 


0 0 0 


STATUS BITS 


WORD 2 
(ODD INPUT CHANNELS) 


Figure 3-3. Data Word Structure 


3. 5 SCALING ALGORITHM 


There are basically two types of numerical values a variable can 
attain during the execution of a computer program. The first type, 
integer values, is generally used where whole numbers are either neces- 
sary or convenient. The second type, real values (commonly called 
floating point numbers), are generally used for all computational values 
where greater accuracy is needed. On most computer systems the limits 
for both real and integer values are sufficiently large to alleviate any 
problem due to overflow; i. e. , to generate a number larger than the 
machine can comprehend. On the CDC 469 mini -computer the limit for 
integer values is sufficiently large, but real numbers must adhere to the 
following constraints: 


C s Real no. < 1. 


Obviously, many computations will not satisfy this constraint, and, there- 
fore, scaling each number becomes a necessary function. 

To ease the burden of scaling, every real number is represented in 
core as two words. The first word represents a floating point number 
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Table 3-2. Data Storage Allocutions 
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satisfying the previous constraint. The second word represents a scale 
factor by which Word 1 is scaled. To obtain the true value for a particu- 
lar variable, the following operations must be performed: 

1) Evaluate Word 1 by real conversion from octal to 

decimal. Note: the decimal place is implicitly placed 
to the right of the left-most digit; the first digit (left 
most) represents either negative or positive numbers 
(0 = positive, 1 = negative). 


Core 

Octal 

Decimal 

Representation 

Value 

Value 

040000 

0.40000 

0. 50000 

060000 

0. 60000 

0.75000 

120000* 

0. 60000 

-0. 75000 

Evaluate Word 2 by integer 
decimal. 

conversion 

from octal to 

Core 

Octal 

Decimal 

Representation 

Value 

Value 

000003 

3 

3 

000010 

10 

8 

177773* 

-5 

-5 

177432* 

-346 

-230 


3) Evaluate the real value of Word 1 and Word 2 by the 
following algorithm: 


X 


(Word 1) x 2 


• (Word 2) 



I 


(See Table 3-3 for further example. ) 




Tables 3-3. Conversion Examples 



Core 

Representation 

Octal 

Value 

Decimal 

Value 

Total 

Value 

Word One 

060000 

0. 60000 

0. 75 

12. 00 

Word Two 

177774 

-4 

-4 


Word One 

130000 

-0. 50000 

-0. 625 

-640. 00 

Word Two 

177766 

-12 

-10 



Negative numbers are represented in core by their "2's" complement. 
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All real values represented by two words can be evaluated in this 
manner. The only exception is the natural logarithm computation. Its 
scaling system is described in the applicable logarithm subroutine 
description, Section 4. 1. 

One of the inherent problems with real number computation is the 
lack of significant figures associated with single precision operations. 
Since single precision offers only five significant figures (a function of 
the CDC 469' s 16- bit word) it becomes necessary to rescale at a fairly 
frequent interval to assure a maximum number of significant figures. 
The scale factor (Word 2) is actually the number of left or right shifts 
performed to obtain this maximum number of significant figures. A 
more complete description of the scaling process can be found in the 
SCALE subroutine description in Section 4. 



SECTION IV 


AUXILIARY SUBROUTINES 


The major portion of the calculations are performed by various 
auxiliary subroutines. The driver's or main program*s sole purpose is 
to control the interaction between these auxiliary subroutines. Each sub- 
routine performs a specific task or calculation based on a specific set of 
inputs. This set of inputs, as well as the output, are generally trans- 
ferred to and from the main driver via the registers. The following sec- 
tion describes the functions of each subroutine, the registers used, and 
shows a simplified block diagram. 

4. 1 INIT (Shown in Figure 4-1) 

The CDC 469 computer is a 20,000g word machine with the majority 
of its core being write-protected (non-volatile). It is desirable to have 
the program residing in this area since it will not be accidentally destroyed 
by either a random write or loss of power. The only section that is not 
write-protected is the scratch pad, words 20g to 20008. This scratch pad, 
as the name implies, is used for temporary storage of diagnostic informa- 
tion and calculations. 

The first page of scratch, words 208 to 3778, have the added capa- 
bility of being addressed directly. This is a very fast method of retriev- 
ing information directly from core. For thi6 reason, any frequently 
used constants are stored here to enhance the speed of the program. 

Since the scratch could be destroyed randomly or by power failure, it is 
not wise to store these constants there permanently. Therefore, any 
constants that are to be frequently used via scratch, are first loaded into 
core in the higher write-protected areas, then at program startup, they 
are transferred to page zero. This is the purpose of INIT. One now 
has the advantage of permanently stored constants in the write-protected 
memory, as well as the added advantage of direct addressing of these 
constants via page zero. 

Calling Sequence : 

SR J INIT, 17 

Registers Used ; 

R2 = Location of data (page zero) 

R3 = Location of data (write-protected memory) 

R4 = Temporary 
R5 - Loop Counter 
R 17 - Exit 

Required Routines ; 

None 
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CORE LOCATIONS 2000 g TO 20000 g ARE WRITE PROTECTED 


Figure 4- 1 . IN1T 


4.2 ZERO (Shown in Figure 4-2) 

The purpose of the routine is to "zero" the scratch pad (location 
20g to 17778)- Since the computer is designed by hardware to start exe- 
cution at location 400g, and the program resides at 2000g f a subroutine 
jump (SRJ) must be re-stored at this location after core has been zeroed. 
This will immediately transfer control to the first word of the program. 


Core Location 
(Octal) 


Contents 

(Octal) 


400 400 

401 2000 


Calling Sequence : 


SRJ ZERO. 17 


Desc ription 

SRJ — Subroutine Jump 
Transfer address 


Registers Used : 

R 1 = 00000 

R2 = Memory scratch address ( 1 7 g to 2000g) 




0 '.‘ 
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ZERO 



Figure 4-2. ZERO 

R3 = 1777 8 (2000g- 1 ) 

R 17 - Exit 

Required Routines : 

None 

4.3 IL1 (Shown in Figure 4-3) 

At 1-second intervals during the execution of the program, a hard- 
ware generated interrupt alerts the computer that data are available on 
the eight input channels. Once the dat* are available it is immediately 
read into memory by this routine. 

The data inputs are simple binary integers. A single 16 bit word 
has a maximum integer value of 77777g. The data inputs may range 
from 0 to l,000,000jo counts. This presents a slight problem since the 
value of the counts received is greater than the capacity of a single word. 
Therefore, two words are used for data inputs. The data format is 
described in Section 3.4. For calculation purposes, this integer value 
for each channel is converted to a simple floating point number, referred 
to as the scaled value. A complete description of the scaling algorithm 
is given in Section 3. 5. The purpose of this routine : s *hen to receive 
the counts from each of the four channels, convert them from two integer 
words to a double precision floating value, and store for further use. 
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Figure 4-3. IL1 


Registers Used : 

R6 = Temporary 
R12 = Input Word 1 
R13 = Input Word 2 
R14 = Temporary 
R15 ■- Masking element 
R16 = Exit 

Calling Sequence : 

Since IL1 is initiated upon a hardware interrupt, nc calling sequence 
is available. 

Required Routines : 

SCALE 

4.4 ZS (Shown in Figure 4-4) 

This routine generates the actual ray travel distance through the 
fuel in centimeters. These calculations are based on the number of 
counts received by the detectors under empty tank conditions and a gamma 
absorption coefficient. Since the standard deviation of the count rale is 
essentially the square root of the counts received, it is conceivable the 








1 


V 


counts at low mass readings will exceed the counts received at empty 
tank calibrations, in this situation the counts received are set equal to 
the calibration inputs. All calculations in this routine are based on the 
scaled inputs as is the evaluation of the logorithm. (See Section 5. ) 


Algorithm : 


Z = 




where 


jjl = absorption coefficient (cm“*) 

N. = counts received 
1 

= counts received with an empty tank 
Z = gamma ray travel distance through the fuel (cm) 
Calling Sequence : 


SRJ 


ZS, 17 

Registers 

Used: 

R2 

- 

Location of data 

R3 

= 

Location of N 

e 

R6 

= 

Counts /answer out 

R7 

- 

Scale factor 

RIO 

= 

N 

e 

Rll 

= 

Scale factor 

R 1 2 


Temporary 

R 15 

= 

Loop counter 

R 17 


Exit 


Required Routines : 

ADDR, DIV, LOG, SCALE 
4. 5 MS (Shown in Figure 4-5) 

This routine performs the actual mass calculations for each of the 
four channels as a function of their respective "Z M length (the gamma ray 
travel distance through the propellant). This computation is actually 
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Figure 4-5. MS 

a polynomial evaluation based on prior data accumulation and polynomial 
curvefit. A representative sample size of mass versus Z length for 
each channel is gathered and fed into an offline computer system with a 
curvefit routine. Each channel's data may be subdivided into two groups. 
For example, a linear curvefit may suffice for Z lengths greater than 
20 centimeters, while a third degree polynomial might be required for 
Z lengths of less than 20 centimeters. A polynomial curvefit is then 
performed for each of these data samples. The result is a set of coeffi- 
cients for each polynomial that represent the curve(s) for each channel. 
These coefficients are then fed into the 469 for the mass calculations 
performed by this routine. 

The polynomial curvefit may attain any degree, although because 
of the word size limitations of the computer, a degree of greater than 4 
would achieve no greater accuracy. Polynomials from one segment or 
channel to the next, need not be of the same degree; the only requirement 
is that the variable NPTS be set to the largest degree of any of the 
polynomials. 
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Calling Sequence : 


SRJ 


MS, 17 

Registers 

Used: 

R2 

= 

Location of "Z" data 

R3 


Location of the coefficients 

R4 

= 

NPTS — number of coefficients 

R5 

= 

Counter 

R6 

= 

Temporary 

R7 

= 

Temporary 

RIO 

= 

Temporary 

R 1 1 

= 

Temporary 

R 12 

= 

Octal mass 

R 13 

= 

Scale factor for mass 

R14 

= 

Z 

R 15 

= 

Scale factor for Z 

R 16 


Exit for scaling routine 

R 17 

— 

Exit 


Required Routines : 

PICK, SCALE 

4. 6 PICK (Shown in Figure 4-6) 

PICK is used exclusively for the purpose of choosing the proper 
curvefit for the mass calculations in subroutine MS. This evaluation is 
based on the variable BKPT, If the calculated ray travel distance, Z, 
is less than BKPT, the first curvefit for the particular channel is used; 
if Z is larger or equal to BKPT, the second curvefit is used. 

Calling Sequence : 

SRJ PICK, 16 

Registers Used : 

R2 = Location of data 
R3 = Location of coefficients 
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PICK 



Required Routines : 
ADDR 


4. 7 FAIL (Shown in Figure 4-7) 

During the execution of the program it is imperative that the four 
channels are generating "feasible" masses. If all four channels are 
operating smoothly and their respective masses are comparatively equal, 
there obviously is no problem. But if a malfunction does arise, this 
routine must decipher which channels are operating properly and what 
course of action is to be taken. The three options listed below are avail- 
able when operating with four source/detectcr pairs: 

1) All channels operating normally 

2) One channel failure — three normal channels 

3) Two or more inoperative channels. 

The system can function properly with three of the four channels 
available with correct data by eliminating the bad data channel from the 
final mass calculations. This is done simply by assigning a weighting 
factor of zero to this particular channel an' 1 a weighting factor of 
0. 3333 to the remaining three channels. (Under normal operating condi- 
u. ' c weighting factor of 0. 25 is assigned to each channel. ) If two or 
more channels are considered inoperable, the system will shut down 
until at least three channels respond to a consistent mass reading. 

The following differences are calculated for each 1- second interval: 


ABS (M 1 - M 2 ) 

(1) 

ABS (Mj - M 3 ) 

(2) 

ABS (Mj - M 4 ) 

(3) 

ABS (M 2 - M 3 ) 

(4) 

ABS (M 2 - M 4 ) 

(5) 

ABS (M 3 - M 4 ) 

(6) 
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Figure 4-7. FAIL 


where: 

ABS = ABSOLUTE VALUE 
KL = Computed mass for channel i. 

These six quantities are then compared to an error criterion, DELM 
(currently 72 grams), and the corresponding equation number is accumu- 
lated in a register if the absolute value of the difference is greater than 
DELM. This accumulation of the equation number then will ide. if'' the 
channel in error. For example, if all six quantities are less than 1 ELM, 
the sum of the equation numbers is zero, signifying all channels are 
operating properly. If channel three is malfunctioning, then Equations 
(2), (4), and (6) will fail the comparison test and the sum of the numbers 
will be 1^10. This sum is then stored into the variable FLAG 1 and passed 
onto subroutine FMASS. One other error indicator is set at this time — 
FLAG 3. This variable represents the actual number of equations that 
failed the comparison criteria. This is also passed onto FMASS. 
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Calling Sequence : 

SRJ FAIL, 17 

Registers Used : 


R2 


Location of 2AM 

R3 

= 

Location of M. 

i 

R4 


Location of M. 

J 

Temporary 

R5 

= 

R6 

= 

M. 

1 

R7 

r 

Scale factor for M. 

RIO 

= 

M. 

i 

R 1 1 

= 

J 

Scale factor for M. 

J 

Temporary 

R 12 

= 

R13 

= 

Temporary 

R 14 


Loop counter 

R15 

= 

Loop counter 

R 16 

= 

Failure Counter 

R17 

r 

Exit 


Required Routines : 

ADDR 

4. 8 FMASS (Shown in Figure 4-8) 

This routine generates the final mass as a function of the four 
channel mass calculations. Depending on the feasibility of the four input 
masses determined by FAIL, this routine weights the masses accord- 
ingly and sums the results to yield a total system mass. The bulk of the 
channel failure selection was done in FAIL and only the two error indi- 
cating variables, FLAG 1 and FLAG 3, are passed onto the routine. 
(FLAG 1 represents the sum of the equation numbers of the failed com- 
parison tests, and FLAG 3 represents the actual number of comparison 
test failures. The following table represents the possible conditions and 
values of FLAG 1 and FLAG 3, 


P 


<*• 


\ 


1 , 
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Figure 4-8. FMASS 



Condition 

FLAG 1* 

FLAG 3* 

No Channel Failures 

0 

0 

Channel 1 Failure 

6 

3 

Channel 2 Failure 

10 

3 

Channel 3 Failure 

12 

3 

Channel 4 Failure 

14 

3 

2 or more Channel Failures 
Octal values 

*15 

>5 


The first function of FMASS is to determine how many channels are 
operating properly. FLAG 3 is used for this purpose. FLAG 3 equal to 
zero indicates all channels are operating properly, and the final mass is 
generated from all four channels by weighting them appropriately 
(wti = 0. 250). FLAG 3 greater than three indicates two or more chan- 
nels are inoperable, and the system cannot generate a mass with suffi- 
cient confidence. The display is immediately blanked, and the program 
recycles itself until at least three of the four channels are functioning 
properly. 

If FLAG 3 has the value of three, a single channel failure is 
detected, and the specific channel must be identified. FLAG 1 is used 
for this purpose, and a minimal amount of comparison logic is needed 
to establish which channel is at fault by the use of the above table. Once 
the malfunctioning channel is identified, it is simply a matter of assign- 
ing the proper weighting factor to each mass to generate a total system 
mass. For example, assume FLAG 1 = 12; referring to the table, it is 
found that channel 3 is not operating properly. Therefore, the weighting 
factor for channel 3 is Bet equal to zero, and the other three weighting 
factors set equal to 0. 3333. Therefore, the total mass is equal to 

+ M2 + M4. From the table, if FLAG 1 = 0, all channels are oper- 
ating correctly, and the respective weighting factors for the four chan- 
nels are set equal to 0. 25. 

The total system mass is a scaled floating point number that is 
now converted to a more meaningful octal integer constant and stored in 
core for later octal to decimal conversion. 

Calling Sequence : 

SR J FMASS, 17 

Registers Used : 

R1 = Location of final mass 

R2 = Location of M. 

1 

R3 = Location of weighting factors (Wt^) 

R4 = Temporary 

hi ) : m/1.,1 ij;l( n y oj,' j 

CUUUiNAL PAbE IS POOR 
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= Counter 


M. * W r l ; 
1 1 


= Scale factor 


Sum of M. * Wtj 


= Scale factor 


Temporary- 
Temporary 
Failure flag 
Temporary 


= Exit for scale 


Required Routines: 


ADDR, SCALE 


4.9 CQNV (Shown in Figure 4-9) 


This is a simple utility routine used to convert the octal integer 
mass to a BCD (Binary Coded Decimal) sing] word output. This can 
then be transfer - directly to the display console indicating an integer 
decimal mass in grams. The algorithm for conversion is as follows: 


ABCD - Octal number to be converted 


WXYZ - BCD word (Binary Coded Decimal) 


WXYZ = A -8 3 + B-8 2 + C-8 1 + D*8° 


Note: For octal representation, three binary digits represent an octal 
digit. For BCD representation, four binary digits represent a 
decimal digit. For example, 


101 = 5 g ; 0101 


BCD 


Similarly: 


001 001 = llg = 9 10 ; 1001 


Q 

'BCD 


Calling Sequence: 


CONV, 17 


Registers Used: 


Location of M. 


R6 - Temporary 
R7 - Temporary 
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Figure 4-9. CONV 


RIO = Number to convert (octal) 

R 1 1 = 00000 

R12 = Conversion factor 

R 1 3 = C0000 

R14 = Output — BCD 

R15 = Temporary 

R16 = Loop counter 

R 17 = Exit 



DOUT is a utility routine used to display the four count readings and 
the mass to the teletype in octal notation. Since the system teletype is 
only capable of printing a line every 2. 5 seconds, the routine should not 
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be called during the operating portion of the gauging system. It should 
rather be used for tape playback only. The variable GOTEL (Location 
20258) governs this option as below: 


Core Location 

2025g 

2025 g 

Algorithm: 


Contents 


Comment 


014000 Enable: Teletype Output 

014001 Bypasses Teletype Output 


1) Input character (from Storage) 

2) Waits for Teletype not busy 

3) Print Character 

4) Repeat sequence (number of iterations 544). 


Calling Sequence : 

SR J DOUT, 17 


Registers Used : 

R1 = Location of data 

R2 = Loop counter (4 channels) 

R3 Loop counter (7 characters /channel) 

R4 = Masking element 
R5 = Temporary storage 
R6 Data register 
R7 Data register 
R 1 3 - Zero — ASC code (260B ) 

R14 - Carriage return (212B) 

R15 * Line feed (215B) 

R 1 6 - Space (240B) 

R17 = Exit register 

4. 1 1 LOG (Shown in Figure 4-11) 

LOG is a s.ngle precision routine written to evaluate the natural 
logarithm of a sealed floating point input. The input, previously scaled 
(or rescaled) by routine SCALE, consists of two words. Word 1 repre- 
sents a floating point number between 0 and 1. Word 2 represents a scale 
factor (a power of 2). Therefore, X, the floating point input, can be 
represented by: 


X = W • 


,K 


(4-1) 
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Figure 4-11. LOG 
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where: 

W = WORD 1 (MANTISSA) and 0 < W <1 
K = WORD 2 (SCALE FACTOR) and -®< K <® 

Since 

In (X) = In (W) + K * In (2) 

and noting the series expansion for In (1 + Z): 

a , 4V n -l n 

E l- 1 ) Z 

n 

n=l 

for 

| 2 | <• 

Z t -1 

and making the substitution, W = 1 +Z, Equation (4-3) becomes 


In (W) = £ 
n= 1 


« (-l) n_1 (W-l) n 


(4-2) 


(4-3) 


n 


(4-4) 


for 


0 < W <2 

Equation (4-4) is a valid expression for In (W) except at the lower 
limit where W = 0. However, a simple "zero-check" of the channel inputs 
causes the series evaluation to be by-passed in the event W = 0 since this 
condition indicates that counts are not being received (i. e. , the particular 
channel is not operating). 

Evaluating Equation (4-4) for a finite number of terms, j, yields: 


for 


(. l )"- 1 z n 

In (W) =^|P +K* 0. 69315 +6 (4-5) 

n- 1 


0 < W <2 
Z = W-l 
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and where: 

In (2) = 0, 69315 

The quantity, €, in Equation (4-5) represents the error inherent in 
In (W) due to only a finite number of terms, j, in the expansion. Obviously, 
it is desired to keep the value of € small, but consistent with the overall 
system accuracy. In the nucleonic gauging system, a value of € <2“^ 
was consistently chosen since it corresponds to the last (15th) bit of a 
469 word. In this manner, any additional number of terms in the expan- 
sion beyond the minimum number required to yield € <2"*^ would not 
affect the final result. 

Algorithm : 


let 



1 * 


I 

1 


X = W * 2 K 


0 < W <1 (WORD 1) 


- co < K <® (WORD 2) 


then 


In (X) = In (W) + K * In (2) 


noting the series expansion for In (1 + Z): 


f(.i) n - 1 z n 

n 

n=l 


if |Z|<1 

zT-i 


let 

W = 1 + Z 

then 

W -1 = Z 

substituting yields: 
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M)"' 1 (W-l) n 

ta « W '=2- n 

n= 1 


-1 < W-l <1 or 0 < W <2 


Ln W is found to be a legitimate expression except at the lower limit, 
i. e. , W = 0; at this limit the series is not evalua'ed since this would imply 
counts are not being received and that that particular channel is not oper- 
ating. A simple zero check on the inputs bypasses the series evaluation 
in this case. 

Evaluating ln (W) for n = 20 iterations yields 

In (W) = Z (i + Z (-1/2 + Z (1/3 + . . . 

+ Z (1/17 + Z (-1/18 X Z (1/19 + Z(-i/20) ). . . ) 


now 


ln (2) = 0. 69315 


then 


K * In 2 = K * 0. 69315 
therefore: 


and 


20 

ln (W) = ^ 
n=i 


, 4V n-l n 
(-1) Z 


n 


+ K * 0. 69315 + € 


(0 < W <2 
7 . = W -1) 



n=2l 


(- 1 ) 


n-1 


n 

Z 


n 



Calling Sequence : 

SRJ LOG, 17 
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Registers Used : 

R4 = Dividend (1.0) 

R5 = Loop counter 

R6 = W Word 1 input 

R7 = K Word 2 input scalefactor 

RIO = LOG W Word 1 output 

Rll = K Word 2 output 

R12 = W-l Word 1 (double precision) 

R13 = W-l Word 2 (double precision) 

R14 = -Temporary storage 
R15 = Temporary storage 
R17 = Exit register 

Required Routines : 

None 

4. 12 EXP (Shown in Figure 4-12) 

EXP is a single precision routine written to evaluate the exponential 
of a scaled, two-word input. This routine was written to work in conjunc- 
tion with the LOG routine. The output of LOG is compatible with the input 
of EXP. The logarithm input is assumed to be input in the following format: 

lnX - InW + k In 2 


where 


| 


i 

* 


InW = Word 1 input 
K = Word 2 input 
therefore: 

Y InW ... Kin 2 w 
X = e e = W v 2 

Hence, W must be evaluated by the exponential of Word 1 input and then 
the proper number of right shifts (multiply by 2) performed to obtain the 
final result. The exponential is found by evaluating a series (see algo- 
rithm this section). Sufficient accuracy was obtained by evaluating 
seven terms of the series. 





Figure 4- 12. EXP 
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Algorithm: 


Q3 r 

L£ 


or in our case 


(InW) 


. i + (lnW) 


Evaluating seven terms and expressing the results a little differently 
yields: 


w = e lnW = 1 + InW |l + InW |jr + InW ( jr + • • • 

(£ + 1 " w (tt) ) ”•) 


+ InW hr + InW 


Hence: 


X - W * 2 1 


Calling Sequence: 


EXP, 17 


Registers Used: 


R4 = N (7 iterations) 

R5 = Loop counter 
R6 = InW input Word 1 
R7 = Input scale factor (2^) 

RIO - Word 1 output (double precision) 
Rll = Word 2 output (double precision) 
R 12 = 1/n! 
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R14 = Temporary storage 
R15 = Temporary storage 
R17 = Exit register 

Required Routines : 

None 

4. 13 ADDR (Shown in Figure 4-13) 

ADDR is a single precision routine written to evaluate the sum of 
two scaled input. Before the addition of any two scaled data inputs can 
be evaluated, they both must have the same scale factors. This is neces- 
sary since Word 2 of a scaled variable is actually a power of two by which 
Word 1 is scaled. Similar to polynomial addition, only like terms of the 
same power can be summed. 

The input with the largest scaled factor is adjusted to the same 
scale factor of the other. The adjustments consist of performing the 
required number of right shifts to the smaller of the two inputs (i. e. , the 
largest scale factor). Once both inputs have the same scale factor, the 
addition can be performed. 

Algorithm : 

Assume inputs: Aj, A 2 (WORD 1) 

Bj, B 2 (WORD 2) 


where 


X, = A, « 2 
X 2 = A 2 * 2 


■B 


-B- 


If Bj = B^, then 


Xj + X 2 = (Aj + A 2 ) * 2' Bl 


If B j < B 2 , then 

Xj - Aj 2" Bl 

X 2 " [ A 2 * 2_(B2 " Bl> ] * 2 " Bl = A 3 :: 2 * Dl 

where 

A 3 = A 2 ^ 2 -(B 2 ' B l ) 

But since B 2 - Bj >0, then 

A 3 = A 2 /2 Bz - B 1 
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Figure 4-13. ADDR 

Now one right shift of an octal number is equivalent to dividing by 2; 
then A 3 is equal to A2 shifted to the right (B2 - Bi) times and thus 

X t + X 2 = (Aj + A 3 ) * 2' B 1 
If B 2 < Bj, create a new variable such that 

V =A z 

V - A i 

B r = b 2 

V = B 1 

Now Bj f < and the previous algorithm can be applied: 
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Example: 

Aj = 23164 Bj = 00004 

A 2 = 16732 B 2 = 000C2 

X t = 23164 * 2" 4 
X 2 = 16732 * 2' 2 

Performing two right shifts in Xj, 

Xj = 04635 * 2" 2 

and 


X { + X 2 = (04635 + 16732) * 2" 2 
= (23567) * 2' 2 

NOTE: For illustrative purposes the sign bit is omitted and assumed 

equal to zero (positive). 

Calling Sequence : 

SR J ADDP , 17 

Registers Used : 

R6 ~ First input Word 1 

R7 = First input Word 2 (scale factor) 

R10 = Second input Word 1/Result Word 1 
Rll = Second input Word 2/Result Word 2 
R12 = Temporary storage 

Required Routines : 

None 

4. 14 DIV (Shown in Figure 4-14) 

DIV is a single precision division routine implemented for the 
purpose of dividing a two-word scaled data input by another two-word 
scaled data input. 

Word 1 of both the divisor and dividend represent, as in all scaled 
data, a real number 0 < X < 1; Word 2 represents the respective scale 
factor. Hence, Word 1 of the dividend can be divided by Word 1 of 
the divisor. Since the scale factors represent powers of 2, the quotient 
is evaluated by subtracting Word 2 of the divisor from Word 2 of the 
dividend. 
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RIO = Divisor Word 1 
Rll = Divisor Word 2 
R12 = Temporary storage 
R17 = Exit register 

Required Routines : 

SCALE 

4. 15 SCALE (Shown in Figure 4-15) 

Since the CDC 469 mini-computer offers only five significant octal 
figures (a function of its 16-bit word), after many operations the results 
must be rescaled to retain maximum accuracy. This is accomplished by 
left shifting the result and dropping non- significant zero bits until a non- 
zero bit appears. Since real values are necessarily <1, only left shifts 
need be performed. The number of left shifts is recorded and becomes 
either the new scaling factor or a method of updating the old one. 

This routine serves both an initial scaling function as well as a 
rescaling function. Initially, data are fed into the central processor as 
raw count data. External peripherals accumulate data as an integer result 
and feed this result into the processor as two words. This input must be 
converted to a scaled two-word data format. The first word represents 
a floating-point real number and the second represents a scale factor. 

This is accomplished by the process described above. Finally, the first 
five most significant bits are retained and stored into Word 1 and the 
number of left shifts is stored into Word 2. 

The rescaling procedure is similar to the initial scaling sequence 
in its technique, but differs in two areas. First, the two-word input 
represents a previously scaled number. Word 1 being a real number 
and Word 2 a scale factor. Secondly, the mber of left shifts is added 
to Word 2. 

Both negative and positive numbers can be scaled or rescaled in 
this manner. 

Example: 

let X = A * 2' B 
where A = WORD 1 input = 03561 

3 = WORD 2 input = 0 
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Figure 4-15. SCALE 


The following table represents the values of A and B after each 
left shift. 


Number of Left Shifts 


Sign 

Bit 


0 

1 

2 

3 

4 

5 


0 
0 
0 
0 
0 
1 > 


03561 

07342 

16704 

25610 

53420 

27040 


B 


0 

1 

2 

3 

4 

5 


Denotes negative number. 

Thus, the required scaling is four left shifts (B -4) for which A = 5342 



Calling Sequence : 


• Initial Scaling: 

SRJ SCALE, 16 


• Rescaling: 

SRJ SCALE + 2, 16 
Registers Used : 

R12 = Word 1 INPUT/SCALE FACTOR word output 

R13 = Word 2 or scale FACTOR input/scale FACTOR output 

R15 = Left shift counter 

R16 = Exit register 

4. 16 SORT (Shown in Figure 4- 16) 

SORT is a single precision routine to generate the square root of a 
real input. The technique used is a recursive formula based on the 
rationale that if Yj is an approximation to Y = \ r x, then 



is a better approximation. Various simulation runs hav'e shewn that two 
iterations are sufficient for determining the square root of a 16-bit word 
(error >Z~ 15 ). 

Since Word 2 is the scaling factor for Word 1 and represents a 
power of 2, it must be adjusted by simple dividing by two, or equivalently 
performing a right shift by one bit. To insure an integer result* Word 2 
should be input as an even integer. The adjustment, if Word 2 is odd, 
must also be made to Word 1 to preserve its initial value. 


Algorithm: (Newton Iteration Technique) 

let X = A * 2" B 
where A = WORD 1 input 
B = WORD 2 input 

then 

x I/2_. a 1/2. 2 -B/2 
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B is first checked for an odd or even value to insure that B/2 is a 
whole number. In the event that B is odd, B is increased by unity and A 
is left-shifted one bit to preserve the initial value of X. If B is an even 
integer, no adjustment is required. 

1 /2 

The first approximation for A is the straight line Y = 0. 6866 X 
+ 0. 3443; this is a good approximation over the interval 0. 3 < X < 1. 0 and 
is particularly suited to the floating point arithmetic of the 469 computer; 
i. e. , X is assuredly <1. 0 and, if X is properly scaled, the most signifi- 
cant bit(s) are available and X >0. 4. (Conversely, if the first bit is zero, 
then X <0. 4 and rescaling is indicated. ) Thus, with 

Y q = 0. 6866 X + 0. 3433 

the first two iterations of the Newton technique yield 

= 1/2 (Y q + A/Y 0 ) 

A 1/2 % Y 2 = 1/2 (Yj + A/Yj) 

Calling Sequence : 

SRJ SQRT, 17 

Registers Used : 

R0 = "P" Register 
R6 = Temporary storage 
R7 = Y (I) 

R10 = X Input saved 

Rll = Masking element 

R12 = X Input/answer output 

R13 = Scale factor input/scale factor output 

R17 = Exit register 

4. 17 MTGO (Shown in Figure 4-17) 

This routine is not part of the Zero-G programming task but is 
included both in the documentation and listing because of its data storage 
capability. MTGO performs the interface with a PERTEC 7000 mag- 
netic tape unit. Logic has been incorporated into the program to auto- 
matically write on the magnetic tape the four channel accumulated count 
data along with the calculated mass. The format is the same as is printed 
on the teletype, and all data are octal. Computer channels four (4) and 
five (5) have been implemented for magnetic tape control, and the vari- 
able GOMT is the flag for writing or bypassing the write routines. The 
code generated now inters a zero into GOMT to perform the data storage 
as the program is running. Data are written at every 1-second interval. 
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MTGO 



1. INPUT IS TERMINATED BY 
TYPING "$" 

Figure 4-17. MTGO 


If the magnetic tape is no longer needed, the instruction at location 2023g 
(MTGO) should be changed as follows: 


Core Location Contents 


2023, 


2023, 


014000 


014001 


Comment 

Enables magnetic write 
Bypasses magnetic tape routines 


These routines have been incorporated into the program logic to 
facilitate data reduction for the mass versus "Z" length curvefits and 
can be eliminated at any time. If the mag routines are used and a print 
of the data is needed, only two operations are needed to read the tape. 
Register zero should be set equal to 6221g (MTRD) and register one set 
to the number of the first file to be read (octal). An end of file is written 
automatically each time the program is started from 2000g. 
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The following is a list of the routines and locations that can be 
eliminated if so desired: 


Routine Name 



Support 

Routines 


Core Locations (Octal) 

6000 - 6106 
6107 - 6156 
6157 - 6?.20 
6221 - 6353 


6354 - 6477 


4. 18 PAUSE (Shown i-> Figure 4-18) 


This routine will cause a pause in the program execution. This 
pause is susceptible to all three interrupt levels as well as the DX mode. 
Therefore, if the pause is to be independent of breaks, the interrupts 
must be locked out before entrance to this routine. The program will 
remain in this status until a bit is entered into the console's data register. 
The data register is initially cleared when the routine is called; therefore, 
any bit entered will terminate the pause and proceed with the program 
execution. 

Calling Sequence : 

SRJ PAUSE. 17 

Registers Used : 

R10 = Console input 
R17 = Exit register 

4. 19 STOPR (Shown in Figure 4-19) 


During the execution of the program the data register on the console 
is checked after each mass calculation sequence for a non-zero input. If 
the input remains zero, the program continues generating new mass data. 
If, on the other hand, input is non- zero, the program discontinues mass 
calculations, writes an end-of-file on the mag-tape and pauses for further 
instructions from the console. At this point, two things may occur: either 
total program termination may result from a master clear; or the data bit 




P 
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Figure 4- 18. PAUSE 


previously entered to halt mass calculation may be cleared, signifying the 
following sequence of events: 

• Rewind tape to beginning of last file. 

• Display on the teletype the contents of the last data file. 

• Transfer control to the start of the original program. 

At this time the program has been re-initialized and proceeds to 
print "ready 11 on the teletype and awdits description data input from the 
teletype. The program is now in the normal executing mode. 


Calling Sequence : 

UJ STOPR 



ZERO OUT 

PROGRAMMER'S 

CONSOLE 


CHECK INPUT 
FROM CONSOLE 




Figure 4- 19. STOPR 

Registers Used ; 

R1 = Argument register for MTREAD 
Rll = Console inputs 

4. 20 MESG (Shown in Figure 4-20) 

As the mnemonic implies the routine generates a message to the 
teletype. This message may be stored at any location in core and may 
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Figure 4-20. MESG 


be entered either at compilation time or manually via the console, 
form of input is an 8-bit ASC II code with two characters per word 
typical example is listed below: 


17000 

152325 

TU 

17001 

151316 

RN 

17002 

120317 

- O 

'7003 

143306 

FF 

17004 

120314 

- L 

17005 

147701 

OA 

17006 

142240 

D- 

17007 

147304 

ND 

17010 

151317 

RO 


Turn off Load NDRO 






The 8-bit ASC II code can be found in Table 5-1. The last word 
should be composed of all zeros to flag end of message. 

Calling Sequence : 

SRJ MESG, 17 

CON MSG1 (location of first word of message — 

17000 g ) 

Registers Used : 

R6 = Character string 
R7 = Teletype status 
RIO = ADDR of message 
Rll = Loop counter 
R17 = Exit 

4. 21 ZCHK (Shown in Figure 4-21) 

This routine provides a simple check on the system to verify counts 
from the detectors are being received. If all four of the input channels 
read zero counts, the display is blocked; the total mass is set to zero; and 
the status flag (word 200) is appropriately set to show system failure. 

A more sophisticated check iB performed in FAIL. If any of the channels 
shows a non-zero input, the routine will exit and the program proceeds 
through its normal sequence of calls to generate the total mass. 

Calling Sequence ; 

SRJ ZCHK, 17 

Registers Used : 

R1 = Location of data inputs 
RIO -- Data 

Rll = Zero (for comparison) 

R12 = Counter 
R 17 = Exit 

4. 22 NES (Shown in Figure 4-22) 

The gauging system, like any fine precision instrument, should be 
calibrated periodically. This is necessary for various reasons: (1) the 
source strength varies over a long period of time due to the natural radio- 
active decay; (2) the "alignment" of the source to the detector may change; 
(3) any object (insulation, supports, etc. ) placed between the source and 
the detector will change the counts received. The recalibration is quite 
simple. With an empty tank, the gauging program is run for a few sec- 
onds (just long enough to gather a few data points for a meaningful sample/. 
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Figure 4-22. NES 


The program is then stopped, the tape rewound and the data displayed on 
the teletype. The average of each column is the new N e for the respec- 
tive channel 1, 2, 3, or 4. 

The program can then be restarted and the system will ask if new 
N e 's are to be entered. At this time the average counts just computed 
are input and the system is now recalibrated. 


4 - 













I 




Calling Sequence : 

SRJ NES, 17 

Registers Used : 


R 1 

- 

Location of message 

R4 

r 

Counter 

R5 

= 

Counter 

R6 

= 

Input 

R7 

= 

Temporary 

RIO 

- 

Teletype Status 

R 1 1 

- 

Counter 

R15 

= 

331B (Y - Yes) 

R 1 7 


Exit 


* 


I 


I 


* 


4. 23 IDENT (Shown in Figure 4-23) 

To distinguish one run from another, a software-generated num- 
erical identifier is written directly onto the mag-tape as well as displayed 
on the terminal, at the start of each run. This was found necessary since 
all the runs aie essentially the same format. This identifier automatically 
increments after each run and is encompassed by the '$' symbol for easy 
recognition on the printouts. The identifier is reset to zero during the 
initial program startup. 

Calling Sequence : 

SRJ IDENT, 12 

Registers Used : 

R3 = Location of Counter (NRUN) 

R4 = Counter 

R5 = Temporary variable 

R6 = NRUN 

R7 = MASK 

RIO = 260B - ASC (Zero) 

R 1 2 = Exit j 

v 

R 1 5 = 244B ($) I 




C EXIT ) 


1 XX * N RUN 

Figure 4-23. IDENT 

4. 24 PAUSR (Shown in Figure 4-24) 

With any real time interactive system, such as the one for the gauge 
program, there exists the need to pause the system's computer. This is 
necessary since the computer is much faster than the response time neces- 
sary for the teletype to react to a command. The general use of PAUSR 
is to allow time for the teletype to perform a complete carriage return be- 
fore the start of each message to be displayed. Therefore, this routine 
performs a number of CRS 20, 17 (circular right shift of R17, 20 bits) 
equivalent to 0. 8 second. The CRS 20 instruction is essentially a NOP 
(NO-Operation), since to shift R17 20 bits is tc leave it unchanged. 


Calling Sequence : 

SRJ PAUSR, 17 


REPRODUCIBILITY op the 
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SECTION V 
DATA REDUCTION 


The Zero-G gauging system is comprised of three basic parts: 

(1) the source/detector pairs, (2) the computer, and (3) the interfacing 
hardware. The computer section can be further subdivided into three 
section^: data input, mass calculation, and cu.ta storage/output. The 
computer is controlled by an internal clock that issues interrupts at reg- 
ular intervals of 1 second. Therefore, the er c ire cycle of data input, 
mass calculation, and data storage/output mu o t be completed in less than 
1 second. For this reason the only visible o^cput during each cycle is the 
LED disp’ y of the final mass. Since output to a terminal is too time- 
consuming, all the necessary information for later data reduc tion is writ- 
ten to a mag-tape. Data input and the mass calculation have been dis- 
cussed in the previous sections. It is the process of writing the information 
to the tape and the data reduction which is described in the following pages. 

5. 1 TAPE GENERATION 

For purposes of data reduction and analysis, information collected 
and generated during each 1- second interval must be stored. P?oer tape 
and/or hard copy generations are impossible because of the time limita- 
tions. Since disk options were not available, the only feasible means of 
data storage was writing to mag-tape. This tape generation offers a 
number of different methods of data reduction. 

First, the tape can be stored for any length of time, remounted on 
the system, and printed on the local system teletype. Obviously, the tape 
can also be rewound at the conclusion of any flight and its contents dis- 
played. This enables a quick visual inspection of the iatei for consistency 
as well as system status. 

Secondly, the tape can he physically transferred to a large scale 
computer for major data reduction and evaluation. This option was found 
to be the most convenient and efficient "a/ to do the proper analysis. 

During each cycle of calculations, six values are written to mag- 
tape — the tour channels of count information, a final mass, and a dummy 
variable (not used at this time). All six values are octal digits represented 
n standard 8- bit ASC code. (See Table 5-1. ) It is the ASC representation 
that is actually written to the tape. Data are stored in this manner rather 
*han octal or BCD for compatibility with the system teletype. 

5. 2 TAPE CONVERSION 

The gauging system utilizes a Control Data 469 mini- computer for 
all logic and computations. This is a 16-bit machine with a core capacity 
of 8K. For data reduction, a lar^e-scale Control Data 6500 is used. This 
is a 60-bit machine with a core capacity of 350K. '’’he obvious difference 
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is the word length, 16 versus 60 bits. A program, hereafter referred to 
as REDUCER, was written to convert the 16- bit word that was written to 
the mag-tape to the 60- bit word size in the larger computer. 

Each tape generally consists of one day's flight, and each flight is 
composed of a number of maneuvers. It is during these maneuvers that the 
gauging system is activated, and data are generated. Each file on the tape 
corresponds to a particular maneuver. Each file is subdivided into a num- 
ber of records. Each record represents one interval (1 second) of infor- 
mation. (See Figure 5-1.) 

The first records of every file are dummy records containing a 
maneuver description in alphanumeric form. These records are used 
only for the system teletype and can be ignored by the large data reduc- 
tion machine. Each subsequent rt cord contains 9-60 bit words or equiv- 
tly 56y — 16 bit words, 3-1/3- 16- bit words being packed into each 
oit word. Each 16- bit word from the 469 is packed into one 18- bit 
segment of the larger 60- bit word as it is read from the tape. This is 
characteristic of data transfer between different word size machines since 
the transfer must be made for each octal number consisting of 3 bus. 
Therefore, 16 bits must be transferred to 18 bits, an even multiple of 
three. Zeros are placed in the remaining 2 bits. Terminating each record 
is an ASC code representing a teletype rub-out (377g). These two condi- 
tions, the record format and the rub-out terminator, enable the appropriate 
unpacking. 

As one ASC coded word (18- bits) is unpac ked, it is decoded to one 
octal character. This character may either be a digit or a blank. If it 
is a digit, it is merged with any previous digit encountered after the last 
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TAPE FILE I 



Figure 5-1. File Structure 



blank. Therefore, all the digits between two blanks represent one piece 
o' information. For example, the first sequence of six digits represent 
the counts received from channel 1 (as shown in Figure 5-2). At the 
completion of this unpacking and decoding sequence for this particular 
record, the six octal numbers are now as they were generated internally 
to the 469 mini- computer and are ready to be stored by the large machine. 
This entire process of unpacking and decoding has a net effect of converting 
the data written on the tape to the original format used by the 469, which 
is now compatible with the larger CDC machine. 

This information is then written to a local file, denoted TAPE6, 
which is ultimately used as the basis for any data analysis to be performed. 
Figure 5-2 represents a dump of the beginning of a typical file (parabola), 
and Figure 5-3 represents a listing of the corresponding local file, i. e. , 
TAPE6. 

Figure 5-3 (TAPE6) represents more than six columns of numbers. 
The first six columns are the actual count data written to the tape. The 
next four columns represent the Z-length conversion corresponding to the 
count data for each channel per the following equation: 



where 

Z i = the ray travel distance through the intervening mass (cm) 

p = mass absorption coefficient (cm ) 

N. = number of counts received by channel . 
i 1 i 

N = number of counts received in the absence of an intervening 
mass 

Also note that the first six columns are octal while th : st four are decimal. 
5.3 REDUCTION 

Once the data have been converted and stored on a file, virtually 
any analysis can be performed. A series of programs were written to 
employ different techniques of data evaluation. They are listed on the 
following pages, with a short description for each program. 

Any references made to TAPE6 or TAPE7 are assumed to be the 
output files generated by REDUCER. TAPE6 is the total unedited data 
file, i.e., the actual conversion of the flight data. TAPE7 is a summary 
of TAPE6 and contains only the average values of the count data over each 
parabola and the number of 1 -second gauging intervals (PTS) during each 
parabola (or 1-g gauging interval). See Figure 5-4. 
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Figure 5-3. Sample Output Format of TAPE6 
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5. 3. 1 REDUCER 


This program converts the tape format of the actual flight data tape 
(generated by the CDC 469) to the tape format required by the CDC 6000 
series computer. Two local files — denoted TAPE6 and TAPE7 — contain 
the resulting output of the conversion (as shown in Figures 5-3 and 5-4). 

5. 3. 2 BUPP 

i 

This program operates on the data of TAPE6 (generated by 
REDUCER) and calculates the average values and standard error of esti- 
mates for the Z-length for a running sequence of gauging data throughout 
each parabola. For example, two integers, Nj and N2, are requested at 
the start of execution. The first integer, Nj, defines at which 1- second 
gauging interval data are to be considered, whereas the second integer, 
N 2 , defines at which 1-second gauging interval beyond which data are not 
considered. In short, and N£ define a sequence of consecutive 1- 
second gauging intervals upon which the following operations are 
performed: 


Zj 


= Z 2 z.. 

i=N 1 l J 

v~v 


1 


Standard Error of Estimate = 



After the indicated operations, Nj and N 2 are then each advanced by unity 
and the operations repeated to form the "running" sequence until N 2 equals 
the last 1- second gauging interval of that parabola. 


5. 3. 3 BUPP4 


This program operates on the data of TAPE6 and calculates the 
average value and the standard error of estimate for the Z-lungths over 
one fixed consecutive sequence of 1- second gauging intervals. BUPP4 
was used primarily to delete the early and late portions of the parabola 
where the propellant was in a (near) one-g configuration and provided the 
"summary" data used as input to ZERO FIT. 

5. 3. 4 ZEROF1T 


The program which relates the calculated Z-lengths to the actual 
mass of propellant wichin the tank is denoted ZEROFIT. ZEROFIT is 
basically a curve- fitting program which utilizes an orthogonal polynomial 
fit by a least squares approximation; i. e. , the error (M ac t ua l * M curve )2 
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L 


is minimized by successive approximations to 

M = A + A.Z + A Z 2 + * ' ‘ * 
o 1 2 

The coefficients, A , so determined from ZEROFIT are ther stored 
in the CDC 469 so that propellant mass can then be calculated in real time 
in the flight environment. 
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SECTION VI 
SYSTEM COMPONENTS 

The following major components are required for the development 
and implementation of the software for the Zero-G operating system. 
(Figure 5-1 presents a functional relationship between these components. ) 

• Control Data Corporation 469 (I) Micro-computer 

• Pertec Magnetic Tape Unit 7X20 


a) 7 track 

b) 7-inch reel 

c) 25 ips 

d) 556 Bpi 


• Control Data Corporation Programmer Console 


ID No. 460-02 Series A 

Part No. 69012100-02 Series 9R 

• TI 733 KSR W/RS232 Interface (or ASR 33 TTY) 


• Teletype Interface 

• Cathode Ray Tube - BALL (MIRATEL) 
6. 1 DIGITAL PROCESSOR 


The Control Data 460 Digital Computer is a fourth-generation, state- 
of-the-art metal oxide semiconductor (MOS)Zlarge scale integrated (LSI) 
digital computer designed to operate reliably in areas of applications 
where small size, light weight, low power, and existence in a severe 
environment are imperative. Physical and electrical characteristics of 
the 469 computer are presented in Table 6-1. 


6. 2 MEMORY 

The memory used in the 469 processor is constructed with plated 
wire elements and is typically referred to as PWM (plated wire memory). 
This type of memory allows random access, is word organized, and does 
not require rewriting after a read cycle. Because of this nondestructive 
readout characteristic, the operating mode of the system is much faster. 
The basic word length is 16 bits with a read cycle time of 1. 6 micro- 
seconds, a write cycle time of 2.4 microseconds, and an access time of 
900 nanoseconds. 
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Table 6-1. CDC 469 Computer Description 


Word Length 

Memory Size 

Interface Levels 

Weight 

Volume 

Power 

Cooling 

Central Processor 


16 bits 
8, 192 words 
0 to 5 volts 
2. 75 pounds 
60 cubic inches 
15 watts 

Conduction and Radiation 


• General Purpose, Binary, Parallel, Single Address 

• Logic: High Level P-MOS/LSI 

• Hardware Registers: 16 

• Arithmetic: Fractional, Fixed Point, Two's Complement 

• Clock: 1 MHz 

• Instructions: 42 Total 


Input /Output 

• Input Channels: ONE, 16 bit, parallel 

• Output Channels: ONE, 16 bit, parallel 

• Interrupt Levels: 3 Plus Direct Execute 

• Serial I/O: 1 IN: 1 OUT 


6. 3 CENTRAL PROCESSOR 

The operational registers for the 469 central processor reside in a 
set of MOS/LSI devices called the register file. This register file con- 
tains sixteen 16-bil registers composed of four MOS/LSI devices, each 
containing eight words of 8 bits each. The file registers are addressable 
as the first 16 words of memory for upward references. However, instruc- 
tions cannot be executed from the register file. 

The central processor of the 469 computer has a basic word length 
(arithmetic register) of 16 bits. For double precision, the word length is 
32 bits (two arithmetic registers). All arithmetic instructions operate in 
a two's complement, fixed point, fractional mode. The instruction rep- 
ertoire is made up of 42 instructions, of which 41 are 16 bits long and 
one instruction, SRJ (subroutine jump), is 32 bits long. This instruction 
repertoire consists of single and double precision add and subtract, s : nglo 
precision fractional multiply and divide, logical operations, load, iti re 
single and double precision shifts, subroutine jump, unconditional jur; '/« , 
interrupt control, input/output (I/O), and functional decisions (condit v’P*l 
jumps and skips). 
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Register assignment is a function of the current interrupt lev -1. At 
a minimum level 0, register "zero" is committed as the program address 
(P) register, and register "one" is selected for the index register. Reg- 
ister zero for the P register is mandatory, although any register can be 
defined as the normal mode index register by using the load index select 
(LIS) instruction. There exist three interrupt levels as v/ell as a mrect 
execute (DX) level. 

Figure 6-1 presents the functional relationship of Control Data 
Corporation 469 computer and the peripherals requisite to software 
development. 

6. 4 PROGRAMMERS CONSOLE 

The programmer's console is connected directly to the central proc- 
essor and displays instantaneously the contents of the register files. The 
user has the option of controlling the execution of the object code either 
under a step-by-step execution mode or a normal execution mode. If a 
certain instruction need be changed, the new instruction and the related 
data word can be entered onto the console and entered directly into the 
object code from the console. 

6.5 TI 733 KSR 

For larger code modifications the 733 KSR can be used. This tele- 
type, similar to the terminals available for TRW/TSS timesharing, gen- 
erally has the function of controlling the entire system. The teletype is 
connected on a controller interface apparatus that provides the necessary 
interface between the more universal teletype unit and the specific central 
processor. 

6.6 MAG-TAPE 

The PERTEC magnetic tape drive serves two purposes. First, via 
the C DC 6600, a compiled version of the program can be loaded directly 
into central memory for later execution by the console. Secondly, once 
the program is operational, the count data generated by each channel and 
the total system mass can be stored on the magnetic tape for each selected 
counting interval. These stored data can then later be batch-processed to 
aid in the data reduction and/or .rve fit purposes for calibration. 
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